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ABSTRACT 



Aims. We present positions and velocities for 796 planetary nebulae (PNe) in the Fornax Brightest Cluster Galaxy NGC 1316 (Fornax 
A). The planetary nebulae and existing kinematics are used to explore the rotation of this merger remnant and constrain dynamical 
models. 

Methods. Using FORS2 on the VLT, the PN velocities were measured using a counter-dispersed slitless-spectroscopy technique that 
produced the largest-to-date sample outside of the Local Group. Spherical, non-rotating, constant-anisotropy Jeans models were con- 
strained by observations of the planetary nebulae and existing integrated light spectra. 

Results. The two-dimensional velocity field indicates dynamically-important rotation that rises in the outer parts, possibly due to the 
outward transfer of angular momentum during the merger. The modeling indicates a high dark matter content, particularly in the outer 
parts, that is consistent with previous estimates from dynamical models, lensing and stellar population models. 
Conclusions. The exceptionally large sample of PN velocities makes it possible to explore the kinematics of NGC 1316 in detail. 
Comparing the results to other early-type galaxies like NGC 1399 and NGC 4594 (M104, Sombrero), NGC 1316 represents a transi- 
tion phase from a major-merger event to a bulge-dominated galaxy like NGC 4594. 

Key words. Galaxies: individual: NGC1316 - Galaxies: elliptical and lenticular - Techniques: spectroscopic 



1. Introduction 

Early-type galaxies hold an important place in the evolution of 
structure in the Universe, but their formation and development 
are not well understood. The halos are of particular interest be- 
cause of the high dark-matter content and the link between the 
orbital structure and the formation history. However, these re- 
gions are difficult to observe because integrated light spectra 
fail to reach sufficient signal-to-noise at large radii. This chal- 
lenge has been addres sed using novel tec hniques on increas- 
ingly large telescopes dProctor et al.l [200')). IFU spectroscopy 
( Weiimans et alj2009[). and discrete tracers like globular clusters 
dSpitler et al.ll2006t ISchuberth et al.ll2010t) and planetary nebu- 
lae. 

Planetary nebulae (PNe) provide a power ful tool for ex 
ploring elliptica l stella r systems beyond 2 - 3R» (|Hui et al.l 
Arnaboldi et al.l Il998b iPeng et al.l 12004 iMendez et alJ 



1995 



2001 



2008; iDohertv et all 120091; ICoccato et alJ l2009t iMcNeil et alJ 



2010). They are bright enough to be detected as individual ob- 



j ects because they emit up to 15% of their ligh t at [OIII] 5007 A 
(iDopita et al.l 119921; ISchonberner et all 120101) . Moreover, their 
spatial distribution follows the stars, so their measured kinemat- 
ics can be co mbined with integrate d-light spectroscopy from the 
inner regions dCoccato et al.ll2.009f) . Using current 8m-class tele- 
scopes, it is possible to generate large enough catalogs of plan- 
etary nebula velocities to constrain dynamical models of nearby 
galaxies at halo radii. 

We present a planetary-nebula survey of the galaxy NGC 
1316. It was selected as a nearby, merger-remnant archetype. 



* Based on observations made with the VLT at Paranal Observatory 
under program IDs 070.B-0089A 



NGC 1316 is a radio source (Fornax A) with some similari- 
ties ( ISchweizeriri980l) to NGC 5128 (Cen A), the subject of 
one of the earliest a nd most successful planetary-nebula sur- 
veys (Hui et al. 1119951). Previous work on NGC 1316 includes a 
survey dArnaboldi et al.ll 1998b which detected 43 PNe. It is the 
brightest cluster galaxy (BCG) in Fornax and the dominant el- 
liptica l in the infalling subcluster described bv lDrinkwater et all 
J200ll) . It is assumed to be the remnant of a merger of two disk 
galaxies based on the i ntermediate age globular cluster system 
dGoudfrooii et al. 2001) and its dusty, unsettled appearance. 

We present PN observations of the galaxy in Sect. [2] Sect. 
[3] covers the creation of the PN catalog including the separa- 
tion of the overlapping NGC 1316 and NGC 1317 populations. 
The results are presented in Sect. |4] and our analysis, discus- 
sion of the implications and future work are are found in Sect. 
[5] Throughout this wor k, we adopt the m easured distance of 
21.5MpcforNGC 1316 dTonrvet al.|[200lh . so 1' = 6.25 kpc. 



2. Planetary nebula observations 

The planetary nebula kinematics were measured using 
counter-dispersed spect roscopy with FORS2 on the VLT 
(lAppenzeller et al ] |1998I) . As shown in Fig.Q] NGC 1316 was 
observed with 7 fields on 3-5 November 2002. The observing 
conditions are summarized in Table [1] 

The FO RS field-of-view is 6'.8x6'.8 or42.5kpc on each side 
at 21.5 Mpc dTonrvet al.ll2001l) . The spatial scale is 0".25 /pixel. 

We used the 1400V grating with a measured mean disper- 
sion of 0.63Apix~'. The measured FWHM of the arclines in the 
calibration frame was 3.8 pixels. 
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Fig. 1. DSS image of NGC 1316 showing the approximate po- 
sitions of the seven fields observed. NGC 1316 sits at the center 
of Field 1 and NGC 1317 is visible in Field 7. The dispersion 
direction is oriented along the major axis of NGC 1316. The 
actual observed area is slightly smaller because the two counter- 
dispersed exposures of each field do not overlap completely. 

Table 1. Observing log for NGC 1316*. 



Field 


P.A. Label 


Exp. Time 


Seeing 






(s) 


(FWHM) 


NGC1316-1 


W140 


4 x 1080 


1.16" 


NGC1316-1 


E320 


4 x 1080 


1.26" 


NGC 13 16-2 


W140 


4 x 1080 


1.01" 


NGC 13 16-2 


E320 


4 x 1080 


1.13" 


NGC 13 16-3 


W140 


4 x 1080 


0.98" 


NGC 13 16-3 


E320 


4 x 1080 


1.16" 


NGC 13 16-4 


W140 


4 x 1080 


0.93" 


NGC 13 16-4 


E320 


4 x 1080 


0.95" 


NGC 13 16-5 


W140 


4 x 1080 


1.10" 


NGC 13 16-5 


E320 


4 x 1080 


1.10" 


NGC 13 16-6 


W140 


4 x 1080 


1.16" 


NGC 13 16-6 


E320 


4 x 1080 


1.11" 


NGC 13 16-7 


W140 


4 x 1080 


0.81" 


NGC 13 16-7 


E320 


4 x 1080 


0.92" 



* Each field is observed at P.A.=140° and 320°, labeled W140 and 
E230 respectively. 



An interference filter, FILT_503_5+86, was used to select 
emission at the wavelength of the redshifted [OIII] line. The fil- 
ter prevents spectra of continuum sources from overwhelming 
the field, and it allows > 50%-transmission detection of PNe in 
the range of -720 to 2280 km s . A plot of the filter bandpass 
and a discussion of how well it is approximate d by a Gaussian is 
available in Sect. 4.4.2 of iMcNeil etaD (120101) . 

In counter-dispersed spectroscopy, the field is imaged with 
a grism in the light path. This causes continuum sources such 
as stars to appear as streaks in the direction of dispersion. 
Emission-line objects such as planetary nebulae appear as dots. 
Another exposure is taken with the position angle (P.A.) rotated 
180°. The distance between the emission line objects detected 



in the P.A.=0 and P.A. = 180 frames is a function of the object's 
velocity. 

The direction of dispersion is along the horizontal axis (x- 
axis). Both the x- and y-axes contain spatial information due to 
the slitless technique. The grism causes an anamorphic distor- 
tion in the focal plane including contraction in the direction of 
dispersion. 

Although the reduction and calibr ation techniques use d for 
this data are detailed in Sect. 2.3 of IMcNeil et al.l d2010h . we 
summarize them briefly here. The geometry of the field of view 
is mapped using a set of three daytime calibration frames based 
on a drilled mask with a grid of holes. The mask is imaged with 
white light in the first calibration frame. In the second and third 
frames, the light from the mask is dispersed using white and 
HeAr lamps respectively. In combination, the three frames allow 
us to measure the local dispersion, the local bandpass (which 
shifts across the field due to the converging beam and the inter- 
ference filter) and the anamorphic distortion. 



3. Cataloging and isolating the planetary nebulae 

The observations were reduc ed and registered usin g similar tech- 
niques to those detailed in IMcNeil et al.l d2010h . Using regis- 
tered images, PNe were detected by blinking the two counter- 
dispersed exposures. A detection is comprised of a pair of [OIII] 
5007 A emission-line sources (one in each exposure) of approxi- 
mately the same brightness and at the same y-position. The sepa- 
ration between the objects must be consistent with the transmis- 
sion limits of the filter. To estimate and improve the complete- 
ness of the catalog, two independent attempts at the detections 
are undertaken a few days or weeks apart. 

The observations for NGC 1316 were done using FORS2 
instead of FORS 1 . This enabled a more effective calibration set- 
up with a drilled mask instea d of the series of calibration slitlets 
used bv lMcNeil et alJ (120 1 Oh . The distortion corre ction and reg- 
istratio n is the same as the one described in §2.3 of lMcNeil et al.l 
d2010l) . but the calibration positions appear less extended in the 
vertical direction, which makes them easier to measure. 

The astrome t ry wa s performed following the steps described 
in lMcNeil et alJ d2010l) . The position of each PNe is taken to be 
the mean of the detections in the two counter-dispersed expo- 
sures. The astrometric solution is based on the measured posi- 
tions of the stars and the co ordinates of the cor responding ob- 
jects in the USNO-B catalog dMonet et al.ll200 3). Our measured 
uncertainty (rms error of the solution) is about 0."3, which in- 
cludes the difficulty of measuring bandpass-shifted continuum 
sources. 

3. 1 . Velocity calibration 

The drilled-mask calibrations used in these observations are 
more accurate than the slitlets used for FORS 1 data. There is no 
reason to suspect that the NGC 1316 PN velocities are incorrect 
in either a relative or an absolute sense. 

To check the agreement between the 7 fields, we matched 
duplicate detections in the areas of overlap. There were 78 du- 
plicated PNe concentrated in fields that were vertically adjacent 
(Fields 1 & 2, Fields 1 & 3, Fields 4 & 5, and Fields 6 & 7). 
In these four cases, the measured shift between the fields was of 
the same order as the error in the shift, so it was not applied. The 
rms of the differences in the pairs was 42 km s . We estimate 
the instrumental error of this sample to be 42/ V2 «30 km s _1 . 
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Fig. 2. Comparison of the veloc ities in this catalog and the sam- 
ple from lArnaboldi et all (1 19981) for the 23 PNe in common. The 
previous measurement has been subtracted from the newer one 
on the vertical axis and the errors have been combined in quadra- 
ture. The median difference is 1.08 km s _1 . The scatter in the 
differences (70.7 km s _1 ) is explained by the combined instru- 
mental errors of the two samples. This indicates that the absolute 
velocity calibration of this sample is accurate. 



We compared the measured PN velocities from thi s sample 
to those previously measured bv lArnaboldi et all (Il998l) . 23 PNe 
are in common between the two samples. The median difference 
between the two measurements is 1 .08km s . The scatter in the 
difference is 70.7 km s _1 , which is roughly expected from the 
combination of the instrumental errors of each sample (30 km 
s _I and 60 km s ). The differences are plotted in Fig. [2] We 
conclude that our absolute velocity calibration is reliable. 

3.2. Decomposing the NGC 1316 and NGC 1317 
populations 

The NGC 1316 PN sample shows localized contamination from 
the nearby Fornax galaxy NGC 1317. It is an SABa galaxy al- 
most directly North of NGC 1316. Their projected separation 
is 377", and their systemic velocities are only 181km s _1 apart 
(NGC 1316 at 1760 km s 1 from Longhetti et al. 1998 and NGC 
1317 at 1941 kms -1 from deVaucouleurs et al. 1991). Due to 
NGC 1316's rotation, the velocity separation near NGC 1317 is 
even less. Although it is a much smaller galaxy, we expect there 
to be contamination from NGC 1317 in the field that covers its 
nucleus, Field 7. The decomposition was applied to Field 7 alone 
because it is the only field significantly influenced by NGC 1317. 
The decomposition of the NGC 1 3 1 6 and NGC 1 3 1 7 PN samples 
follows a similar two- fold membership-ass ignment technique to 
the one described in iMcNeil et"ai1 (l2010h for NGC 1399 and 
NGC 1404. Each PN is assigned a probability of membership 
in NGC 1316 based on its velocity and the surface brightness 
of the galaxy at that position. The probabilities are normalized 
such that all objects in the catalog belong to either NGC 1316 or 
NGC 1317, but membership is assigned only if there is greater 
than 90% confidence. 




1000 1200 1400 1600 1800 2000 2200 2400 



Velocity (km s A -l) 

Fig. 3. Velocity distribution and histogram of the PNe from the 
F7 observations. The velocities are represented in the histogram. 
The smooth function represents the maximum-likelihood fit to 
the distribution — it is parametrized as the sum of two Gaussians 
multiplied by a Gaussian representing the filter. 

3.2.1. Velocity criterion 

The first criterion for membership is based on velocity. Because 
the velocity separation is so small, and NGC 1317's systemic ve- 
locity is so near the edge of the filter, we considered the shape of 
the filter to be significant in shaping the overall velocity distri- 
bution. Using a maximum-likelihood technique, we fit the Field 
7 velocity distribution with a sum of two Gaussians (represent- 
ing the two galaxies) multiplied by a Gaussian (representing the 
filter). The fit to the velocity distribution is shown in Fig. [3] The 
fitted means of the two galaxies are found at 1826.86km s _1 and 
1936.08 km s co mpared to the observed values of 1760 km 
s -1 an d 1941km s -1 ( lLonghetti et alJfT998hlde Vaucouleurs et al.l 
[1991) . Accounting for rotation in the outer parts of NGC 1316, 
we find this fit is consistent with our physical understanding of 
the system. 

We then reverse the technique to calculate the probability of 
membership of each velocity in the NGC 1316 Gaussian mul- 
tiplied by the filter versus the NGC 1317 Gaussian multiplied 
by the filter. The probabilities were normalized such that each 
PNe must either be a member of one galaxy or the other. The 
probability of membership in NGC 1316 defines our velocity 
membership criterion. 

3.2.2. Spatial criterion 

For the spatial criterion, we project our PNe onto the major axes 
of the two galaxies and estimate the light contributions from 
each galaxy at the location of each PN. We can use this selec- 
tion m ethod because PNe trace the stellar light dCoccato et al.l 
120091) . 

We estimate the major-axis radius of each PN assuming a 
fixed ellipticity of e = 0.33 based on isophotes from the DSS 
image of NG C 1316 and e = 0.17 fo r NGC 1317 based on mea- 
surements by IWozniak et alj ( 1 1995b . The position angle of the 
photometric major axis of NGC 1316 is assumed to be constant 
at 50°. For NGC 1317, we estimate PA.=150° based on Fig. 2 
m IWozniak et all d 19951) . but this is only valid because we are 
interested in values outside ~ 16". 

We compared the light contribu tion from each gal axy by us- 
ing the B -band light profiles fr om ISchweizeri (Il980h for NGC 
1316 and IWozniak et al.l d 19951) for NGC 1317. To estimate the 
surface brightness profile beyond the limits of these observa- 
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tions, we fit the outermost 14 observations from Schweizer with 
a Sersic profile and the outermost 15 observations from Wozniak 
et al. with an exponential profile. The probability of membership 
based on position is the ratio of the surface brightness of NGC 
1316 to the surface brightness of NGC 1317. 

3.2.3. Selection of the NGC 1316 sample 

To select the final sample for NGC 1316, we multiplied the ve- 
locity and spatial criteria. The final probabilities are normalized 
so that the total probability of membership for each object is one. 
However, objects are only assigned membership in a galaxy if 
they have more than 90% probability of belonging. This yields 
a catalog of 796 PNe bound to NGC 1316. 69 PNe are members 
of NGC 1317 and 43 are ambiguous. 

3.3. Contamination from other sources 

Based on the two-fold decomposition technique in Sect. 13.21 we 
are able to identify and remove PNe that belong to nearby galax- 
ies. Our observations also include possible contamination from 
other sources. Most background galaxies can be identified as 
non-planetary nebulae because they are associated with contin- 
uum, they are spatially extended or they have double emission 
lines. These contaminants are removed before the sample is as- 
sembled. 

However, some emission-line sources cannot be distin- 
guished from PNe — they are most likely unresolved Lyman-a 
galaxies at z = 3.1. These background objects show a veloc- 
ity distribution centered on the filter instead of the galaxy's sys- 
temic velocity enabling us to estimate how much contamination 
is present in the catalog. In the case of NGC 1316, these objects 
would have lower velocities on average than the PNe from the 
galaxy. However, they are extremely difficult to remove because 
their spatial distribution is approximately flat over the field of 
view. Contamination from these background galaxies is mini- 
mized by avoiding faint objects in the selection process because 
their equivalent widths are low and t hey are less likely to be 
continuum-free dArnaboldi et al.ll2002h . In the case of the NGC 
1316 sample, contamination is limited (see the excess in the his- 
togram at low velocities in Fig. |4j> and we are able to effectively 
analyze our sample with methods that are robust to outliers. 

3.4. PN catalogs 

3.4.1. Finalizing the catalogs 

The velocity histogram of 908 PNe observed in the NGC 1316 
fields is shown in Fig. |4] This is an exceptionally large catalog. 
796 of these objects are bound to NGC 1316 with greater than 
90% probability. The sample also includes 69 members of NGC 
1317's PN population and 43 ambiguous objects. The catalogs 
for NGC 1316 and NGC 1317 are presented in Tables and [3] 
respectively. 

3.4.2. Comparing the efficiency of our detections 

The expected number of detections for our sample is determined 
by the observations and the intrinsic properties of the galaxies. 
We can judge the efficiency of each observing run by scaling the 
size of the samples by these factors. 

This sample is comprised of nearly 5. 5 times as many ob - 
jects as the NGC 1399 catalog presented in lMcNeil et alj(l2010l) . 
Since they are both members of the Fornax cluster, we assume a 




500 1000 1500 2000 2500 3000 
Velocity (km s ') 

Fig. 4. Velocity histogram of the PNe detected in the NGC 1316 
observations. This includes objects associated with NGC 1316 
and NGC 1317. The red Gaussian represents NGC 1316 and the 
black curve indicates the transmission of the filter in the relevant 
velocity range. Some low-velocity objects are skewing the dis- 
tribution, and we see the suggestion of a second peak near the 
systemic velocity of NGC 1317. (1941 km s -1 according to de 
Vaucouleurs et al. 1991) 

Table 2. Catalog of PNe in NGC 1316 



ID 


RA 


Dec 


wavelength 


Vhclio 




12000 


12000 


A 


km s" 1 


NGC1316 1 


50.47341 


-37.21804 


5039.0 


1924 


NGC1316 2 


50.47824 


-37.23466 


5033.4 


1587 


NGC1316 3 


50.47914 


-37.21044 


5038.7 


1901 


NGC1316 4 


50.48403 


-37.22564 


5035.3 


1700 


NGC1316 5 


50.48860 


-37.24886 


5031.5 


1472 


NGC1316 6 


50.48942 


-37.21301 


5036.9 


1798 


NGC1316 7 


50.49080 


-37.24301 


5036.4 


1764 


NGC1316 8 


50.49086 


-37.20056 


5033.6 


1598 


NGC1316 9 


50.49483 


-37.20033 


5035.5 


1713 


NGC 13 16 10 


50.49716 


-37.23789 


5033.2 


1572 


NGC1316 794 


50.86464 


-37.18436 


5038.1 


1868 


NGC1316 795 


50.86579 


-37.18282 


5039.4 


1947 


NGC1316 796 


50.87387 


-37.21063 


5038.2 


1876 



NOTES - Details on the astrometry are in Sect. [3] Velocities, 
measured from the [OIII] 5007 A line, are corrected to heliocen- 
tric. The velocity error is 30 km s -1 . The full catalog is available 
in the associated electronic material. 

common distance and compute the size of a hypothetical sample 
for comparison. 

Beginning with the 146 PNe detected in NGC 1399, we 
add 55 objects based on symmetry arguments to account for 
the 4 of 9 fields that were not observed. We scale by 2.42 
to account for the difference in brightness. For this we used 
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Table 3. Catalog of PNe in NGC 1317 



ID 


RA 


Dec 


wavelength 


I'helio 




J2000 


J2000 


A 


km s -1 


NGC1317 1 


50.65941 


-37.09871 


5039.8 


1970 


NGC1317 2 


50.66192 


-37.09394 


5040.7 


2026 


NGC1317 3 


50.66277 


-37.10502 


5040.3 


1998 


NGC1317 4 


50.66303 


-37.10728 


5041.5 


2072 


NGC1317 5 


50.66341 


-37.10949 


5040.1 


1989 


NGC1317 6 


50.66537 


-37.08525 


5039.6 


1956 


NGC1317 7 


50.66625 


-37.09976 


5041.6 


2081 


NGC1317 8 


50.66925 


-37.10303 


5039.5 


1951 


NGC1317 9 


50.66943 


-37.08724 


5039.1 


1927 


NGC1317 10 


50.66995 


-37.10643 


5039.5 


1954 


NGC1317 67 


50.71118 


-37.09235 


5039.2 


1933 


NGC1317 68 


50.71183 


-37.09738 


5039.3 


1938 


NGC1317 69 


50.72124 


-37.10676 


5039.1 


1926 



NOTES - Details on the astrometry are in Sect. [3] Velocities, 
measured from the [OIII] 5007 A line, are corrected to heliocen- 
tric. The velocity error is 30 km s -1 . The full catalog is available 
in the associated electronic material. 



2MAS S (ISkrutskie et al.ll2006h K ma gnitudes adjusted accord- 
ing to iKormendv & Tremaine d20TTh . Finally, we consider the 
relation from Buzzoni et al relating the PN number den- 

sity to the color of the galaxies. The overall trend implies that 
redder galaxies have a higher specific frequency. Despite being 
reddened by dust, NGC 1316 is bluer than NGC 1399 implying 
an increase in our predicted sample. There is significant scatter 
in the Buzzoni et al. relation, accounted for by variation in the 
lifetimes of PN populations, and a much steeper relation for ex- 
ternal ellipticals than for Local Group galaxies. The adjustment 
for our predicted sample is between 1.2 and 2.5. 

Adjusting for these differences, we can compare our sample 
of 796 PNe in NGC 1 3 16 to a predicted sample between 588 and 
1216 PNe. We conclude that, although NGC 1316 has a remark- 
ably large sample, it is not unexpected. 

4. Results 

4.1. The spatial distribution 

The spatial distribution of the PNe in NGC 1316 is shown in 
Fig. The overall shape is determined by the layout of the 7 
observed fields in Fig. [TJ and we see rotation along the major 
axis (PA.=64°, see Sect.l4~2l. 

The analysis of the PN kinematics and their use as trac- 
ers in dynamical models are depende nt upon the consistency 
of PNe and stars. ICoccato et al.l (|2009) showed this for a sam- 
ple of elliptical galaxies observed with the Planetary Nebulae 
Spectrograph. In Fig. [6] we show a comparison of the stellar 
surface brightness and the PN number density. The PNe were 
divided into annular bins and the logarithmic projected number 
density in each bin was calculated as 

/"PNe = -2.5 log(AfpNeM) + C, 

where Np^ e is the number of PNe in the bin, A is the area of the 
annulus and c is an arbitrary constant selected for the galaxy (but 
remaining constant for all the bins). 
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Fig.5. Spatial distribution of the 908 objects in the NGC 1316 
observations. North is up and East is left. The centers of the 
two galaxies are marked with black crosses. Surrounding NGC 
1317, the solid, black dots represent PNe bound to NGC 1317 
with >90%probability. The colored dots indicate membership 
in NGC 1316. If a PN cannot be assigned membership with 
>90% probability, it is marked as ambiguous (hollow black 
dots). Ignoring the structure from this companion galaxy, this 
figure shows NGC 1316's rotation oriented along PA. =64° (see 
Sect. 14.21 ). The north-east is rotating away from the observer 
while the south-west is moving towards us relative to the mean. 



The first bin is incomplete because of the void in the center 
where PNe are not detectable. The surface brightness profiles 
are in good agreement until ~ 250" where the influence of the 
dust becomes apparent. Dust strongly affects the integrated light 
observations in the outer parts. The NE major axis is the least 
affected, and, therefore, is the best point of comparison for the 
PNe. 

In the region where the PN sample is complete and dust does 
not obscure the B-band integrated light, the PN number density 
shadows the NE major-axis surface density of the stars. We con- 
clude that the PNe and the integrated light are compatible trac- 
ers of the underlying stellar mass. We will use the two tracers 
together to constrain dynamical models of NGC 1316 in 35.21 

4.2. The kinematics 

A phase-space projection of line-of-sight velocity against radius 
is shown in Figj7] There are some high- and low-velocity objects 
that may be contaminants, but they will not affect the robust dis- 
persion measurements discussed in Sect. 15.11 

We quantify the rotation of the sample by modeling the to- 
tal rotation and the rotation in annular bins. SO galaxies usually 
display rotation, and we are interested to learn to what extent the 
flattening is rotationally supported. The fitted rotation for the full 
sample of PNe is shown in Fig|8] 

For the fit, we assume that the galaxy is spherically symmet- 
ric, the amplitude of rotation is only a function of radius and that 
the axis of rotation is in the plane of the sky dCote et al.l 1200 it 
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Fig. 6. Surface density profiles for the PNe and stars in NGC 
1316. Red and black crosses mark the NE and SW major-axis B - 
band surface brightness profiles respectively (ICaon et al. Ifl994h. 
The h ollow squares show the measurements from Schweizerl 
(fl980h . The agreement with the NE major axis between 200" 
and 500 " indicates that the PNe do trace the light, except where 
dust obscures the observations. 
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Fig. 7. Line-of-sight velocity against projected radius. The 
sample is centered o n the systemic velocity of 1760km/s 
dLonghetti et alj 119981) . Solid red objects have been assigned 
membership in NGC 1317 with >90% probability. The hol- 
low red dots are ambiguous- they could not be assigned to ei- 
ther galaxy with more than >90% probability. There are some 
extreme-velocity objects near the center of the filter (779km s _1 ) 
and around 2700 km s _1 that may be background galaxies. 



IWoodlev et alJl2010h . The velocity as a function of position an- 
gle, 0, is parametrized as 

v P l&] = V sys +Acos(0-0 o ), 

where V sys is the systemic velocity, A is the amplitude of rota- 
tion, and 0o is the PA. of the kinematic major axis. We also con- 
sidered adding a sin(30) term, but the simpler model was pre - 
ferred by a Bayesian Information Criterion test dSchwarzj|1978l) . 

Because we suspect a fit including a free systemic velocity to 
be influenced by the contamination f rom the center of the filter, 
we hold y sys constant at 1760 km s" 1 dLonghetti et al.ll998l) . The 
result is a kinematic major axis at P.A.=64±8° and an amplitude 
of rotation of 85±11 kms -1 . The photometric major axis is at 
PA. =50°. 

Using the rotation fitted to the entire sample as an initial 
guess, we divided the PNe into five radial bins. The PA. and 
amplitude of rotation are modelled for each bin, and the results 
are in Fig. [9] and Fig. [TOlrespectivelv. The position angle of the 
kinematic major axis is offset from the photometric major axis 
by 1.6cr, but still consistent (see Fig. |9j). The amplitude of rota- 
tion rises in the outer parts. This could indicate that some PNe 
bound to NGC 1317 with a systemic velocity of 1941 km s _1 
dde Vaucouleurs et al.l 19911) have stayed in our sample. This sce- 
nario is not supported by the distribution of the PNe in Fig.[8]nor 
the consistency of bin 5's position angle with the inner bins in 
Fig-|9] The rise could also be a real description of the kinemat- 
ics of NGC 1316. This phenomenon could be evidence of the 
outward transfer of an gular momentum resulting from a merger 
dOuinn&Zureklll988l) . 
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Fig. 8. The line-of-sight velocities of the NGC 1316 PNe as a 
function of position angle. The grey dots represent the decom- 
posed PNe sample. The black fit represents rotation of a spheri- 
cally symmetric system with the axis of rotation in the plane of 
the sky. It indicates an amplitude of 85km s _1 and a kinematic 
major axis at P. A. =64° (1.11 radians). 



5. Discussion 

5. 1 . Line-of-sight velocity dispersion 

The line-of-sight velocity dispersion profile is a useful tool to 
illustrate the mass distribution and dynamics of elliptical galax- 
ies. For our measurements, we use a robust technique for mea- 
suring the dispersio n of discrete tracers that is described in 
iMcNeil et all d2010h . Because the calculation is only based on 
objects lying within 2<x of the mean, outliers do not have a strong 
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Fig. 9. Fitted values for the position angle of the kinematic ma- 
jor axis of NGC 1316. The connected black points show the 
kinematic major axis's P.A. known from fitting the rotation in 
five radial bins. The P.A. fitted to the entire sample is 64±8° 
(shown as grey dotted line). The dashed grey line shows the pho- 
tometric major axis at 50°. 
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Fig. 10. The rotation curve of the PNe in the outer parts of 
NGC 1316. The black line traces the amplitude of rotation of the 
PNe in annular bins. At small radii, the gray dots show positiv e 
major-axis rotation measurements from iBedregal et ail (120061) . 
and the black dots show major-axis rotation measurements from 
lArnaboldi etaD <ll998l) . 



effect. This method is resilient against any lingering contamina- 
tion from background galaxies. 

In the case of NGC 1316, we measure the dispersion after 
accounting for the rotation. The robust technique was applied to 
the residuals of the rotation fit in each of the five annular bins. 

The random-motion profile of NGC 1316 is shown in Fig. 
m At small radii, the integrated-light measurements from 
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Fig.ll. Line-of-sight velocity dispersion profile for NGC 1316. 
The hollow circles represen t the integrated light measurements 
from lBedregal et ail d2006l) for the major axis. The PN disper- 
sions (hollow stars) are computed using the robust technique 
on the residuals of the NGC 1316 rotation fit. The instrumental 
dispersion is subtracted in quadrature. The results are resilient 
against outliers and describe the random motion in NGC 1316. 
The two tracers interface well at ~150". 
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Fig. 12. V/cr profile of NGC 1316. The fitted rotation of the PNe 
is divided by the robust dispersion measurement shown in Fig. 
[TTIto illustrate the relative dynamical importance of rotation and 
random motions. For NGC 1316, rotation is more important in 
the outer parts tha n it is in the inner parts. lR e = 109"= 1 1.4kpc 
(ICaon etalJli99l . 



Bedregal et alJ J20Q6) show a falling dispersion profile. As we 5.2. Jeans Models 



move into the halo, the planetary nebulae indicate that the veloc- 
ity dispersion continues to fall off until at least ~ 500". 

V/cr is used to quantify the relative importance of rotation 
and random motions in the dynamics of the galaxy. We divide 
the amplitude of rotation for the PNe by the robust dispersion 
shown in Fig. QT| The V/cr profile for the planetary nebulae in 
NGC 1316 is shown in Fig. [12] Rotation is more important in the 
outer parts than it is in the inner parts of NGC 1316, suggesting 
that angular momentu m transferred outwards during the merger 
(IOuinn&Zureklll988l) . 



To better understand the meaning of the kinematics of the 
PN sample in NGC 1316, we construct a simple dynami- 
cal model of the galaxy. We begin with the Jeans equations 
dBinnev & Tremainell 19871) . assuming spherical symmetry: 



1 d(p^.(r)) B cr 2 r (r) 



P 



dr 



+ 2- 



d<S> 
dr ' 



(1) 



where p is the spatial density, cr r is the radial component of the 
velocity dispersion, (3o is the anisotropy of the velocity distribu- 
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tion (J3 = 1 — -J), which is assumed to be independent of radius 
here, and i> is the total potential. The expressions on the left- 
hand side of Eq. Q] can be related to observationally accessible 
quantities. In this case, we solve for the velocity dispersion and 
correct for the projection onto the plane of the sky. The projected 
velocity dispersion of the model is given as 



31)0 



2 r 
W)J R 



R 2 \ j{r)oi(r)r 
1 -y3o-r . dr, 
r 2 ) yfri^W 



(2) 



where I(R) is the observed surface brightness profile and j(r) is 
the light density which is calculated from the mass density, p, 
divided by a constant mass-to-light ratio, F. 

The velocity dispersion of the model, o> depends on the to- 
tal mass of the system. As a first attempt, we assume a model 
where mass traces the light with a constant mass-to-light ratio. 
We compared the model's projected velocity dispersion pr ofile 
to the observations of integrated light (IBedregal et al.ll2006l) and 
the robust dispersion measurements of the PN sample. The in- 
tegrated light positions were adjusted by -0.32kpc to maximize 
the symmetry of the dispersion profil^B Since this non-rotating 
model does not otherwise account for the energy associated with 
rotation, we constrain our model with the observational quantity 
cr, ,. In the case of the integrated light, this is the line-of-sight 
velocity and the velocity dispersion added in quadrature. In the 
case of the PNe, it is the robust dispersion of the binned PNe 
velocities without adjusting for the rotation. The goodness-of-fit 
was quantified with a reduced ^ 2 -style merit function 



x 1 



N ■ 



i y (gjgtW - o-p(Ri)r 

O"enoi(0 



»tl 



(3) 



where N is the number of data points constraining the fit and 
n is the number of free parameters, <x to t is the observed value 
of sigma at a radius, o- p (Rj) is the projected dispersion of 
the model at R, and cr errol is the uncertainty of the observation 
belonging to that radius. 

The best fitting self-consistent model is shown in Fig. [13] We 
assume that the total mass traces the light with a constant mass- 
to-light ratio. This model has a mass-to-light ratio of 4.OM /L 
and an anisotropy of 0.44. The disagreement between the model 
and the observations in the outer parts indicates the necessity of 
a dark matter component, in this case parameterized as a pseudo- 
isothermal distribution such that 



PdmO) 



DM 



2nG(a 2 + r 2 Y 



(4) 



where cr DM is a velocity dispersion representing the total mass 
of the DM halo and a is a core radius. 

Adding in the dark matter component, we explored a grid of 
models varying over the four parameters /?o, F, <x DM , and a. We 
resolved the anisotropy (J3) to 0.01, the mass-to-light ratio (r) 
to O.1M Q /L , the core radius (a) to 1 kpc and the dark matter 
dispersion (<x D m) to 1km s . 

The best value for our merit function corresponded to a 
model with y8 = 0.50, F = 2.8M /L , a = 5kpc, and cr DM = 289 
km s _1 . The best model is shown in Fig. [14] 



1 iBedregal et al.l d2006l) selected the peak of the dispersion profile as 
their position zero-point. Upon closer inspection, the highest dispersion 
value was an anomaly in an otherwise symmetrical profile. We removed 
the anomalous point and fitted the profile with a symmetrical function 
finding a shift of -0.32 kpc relative to their system. The fit shared the 
same peak as the measurements (with the object representing the previ- 
ous zero-point removed). 
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Fig. 13. Self-consistent Jeans model of NGC 1316. Blue p oints 
are integrated-light dispersions from IBedregal et al.l d2006l) and 
the red points show the dispersions of PNe in annular bins. 
The blue curve represents a model with a mass-to-light ratio of 
4.OM O /L and a constant anisotropy of 0.44. The disagreement 
between the model and the observations in the outer part under- 
lines the importance of DM in the dynamics of NGC 1316. 



300 



250 



200 



150 



100 




1.0 5.0 10.0 

Radius (kpc) 



50.0 100.0 



Fig. 14. Dynamical model of NGC 1316 including dark 
matter. The blue po ints are integrated-light dispersions from 
IBedregal et al.l (120061) and the red points are measured from PNe. 
The model (blue curve) has a constant mass-to-light ratio of 2.8 
Mq/Lq and a constant anisotropy of 0.50. The dark matter distri- 
bution is characterized by <r DM = 289 km s and a core radius 
of 5 kpc. This model does not rotate. 



The modeling of NGC 1316 is limited by the simplifications 
in the modeling technique. The outcome would be improved by 
assuming an axisymmetric distribution that is flattened by rota- 
tion. The mass-to-light ratio might well change over the radial 
range of our models. Allowing the anisotropy to change as a 
function of radius will add flexibility to the fitting and create a 
more accurate model of NGC 1316. 

Despite the inherent assumptions in this model, it presents 
a realistic physical description that is consistent with observa- 
tions and theory. Our mass-to-light ratio of 2.8 M /L lands 
at the low end of values measure d for other galaxies (see Fig 
4.25 in lBinnev & Tremaind (120081) ). The dark matter content of 
NGC 1316 shown in Fig. [15] is towards the high end of esti- 



mates of 5-65% within 1R. (11.4kpc = 109". [Caon et al.ll 19941) 



dGerhard et al.ll2Q01t iTreu & Koop mans 2001 ICappellari etakl 
120061: iThomas et al.ll2007t iBarnabe et al.lt2009h . 
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Fig. 15. Dark matter mass fraction of NGC 1316. By consid- 
ering the dark and luminous components our model separately, 
we can assess what fraction of the mass at each radius is DM. 
This plot is made for the model that best fits the kinematic data: 
/?o = 0.50, T = 2.8M /L O , a = 5kpc, and 0" DM = 289 km s" 1 . 

5.3. Implications for the assembly of early-type galaxies 

Evide nce from the globular cluster system dGoudfrooii et al.l 
1200 ll) indicates that the last major merger occured 3 Gyr ago. 
The resulting galaxy is dusty with dynamically important rota- 
tion. Projecting the evoluti on of NGC 1316 forward in time , the 
dust will settle into a disk (Ivan der Kruit & Freemanll201 lL and 
references therein). A system like this could evolve into an ob- 
ject like the Sombrero Galaxy (NGC 4594, M104). 

NGC 1316 provides a lower limit for the time between a 
major merger and a galaxy like NGC 4594 to form. Likewise, 
it constrains the time since NGC 4594 last underwent a major 
merger to be greater than 3 Gyr. In the context of mass assembly 
within clusters, NGC 1316 tells a different story to the Fornax cD 
galaxy, NGC 1399. NGC 1399 is a slow rotating, possibly triax- 
ial, dominant galaxy. The cD envelope implies that accretion is 
currently dominating its mass assembly. In contrast, NGC 1316 
is in a phase where it has recently undergone a major merger 
and would need to be involved in several more to resemble NGC 
1399. As shown in Fig. [l0]and[TT] NGC 1316's rotation stays 
constant in the outer parts, but the random motions decrease. 
Despite dominating its subcluster, NGC 1316 illustrates an ear- 
lier phase in mass assembly within clusters. 
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5.4. Future work 

The ultimate goal of these observations is to constrain the 
outer parts of dynamical models. We have presented spher- 
ical, non-rotating Jeans models. Future work will focus on 
more sophisticated mo dels including NMAGIC modeling (see 
Ide Lorenzi et al.ll2008l) . 

Acknowledgements. We thank A. Bedregal for generously sharing the major- 
axis kinematics of NGC 1316. We thank W. Harris for discussions regarding the 
rotation of NGC 1316. This work has made use of the NAS A/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under contract with the National Aeronautics and Space 
Administration. 



References 

Appenzeller, I., Fricke, K, Fiirtig, W., et al. 1998, The Messenger, 94, 1 
Arnaboldi, M., Aguerri, J. A. L., Napolitano, N. R., et al. 2002, AJ, 123, 760 
Arnaboldi, M., Freeman, K. C, Gerhard, O., et al. 1998, ApJ, 507, 759 



